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COMPLEXITY OF THE XY ANTIFERROMAGNET AT FIXED MAGNETIZATION 


ANDREW M. CHILDS 1 ’ 2 ' 3 , DAVID COSSET 1 ’ 2 ’ 4 , AND ZAK WEBB 2 ’ 5 


Abstract. We prove that approximating the ground energy of the antiferromagnetic XY model on a simple 
graph at fixed magnetization (given as part of the instance specification) is QMA-complete. To show this, 
we strengthen a previous result by establishing QMA-completeness for approximating the ground energy of 
the Bose-Hubbard model on simple graphs. Using a connection between the XY and Bose-Hubbard models 
that we exploited in previous work, this establishes QMA-completeness of the XY model. 


1. Introduction 

Kitaev pioneered the study of quantum constraint satisfaction problems, where the goal is to approximate 
the ground energy of a Hamiltonian [5], Many examples of such ground energy problems are known to be 
complete for the complexity class QMA, a quantum analogue of NP (see, e.g., [2]). 

In this paper, we focus on computational problems defined by graphs. For example, MAX-CUT is a 
classical constraint satisfaction problem defined by a graph. The goal is to find a subset of vertices that 
maximizes the number of edges between that subset and its complement. This can be rephrased as a ground 
energy problem: it is equivalent to minimizing the energy of the Ising antiferromagnet on the graph, where 
there is a bit for every vertex and a constraint penalizing adjacent bits that agree. Equivalently, we may 
consider a qubit at every vertex and a ZZ interaction for every edge. 

To obtain a genuinely quantum constraint satisfaction problem defined by a graph, we consider interaction 
terms with nonzero off-diagonal matrix elements. Natural choices include the antiferromagnetic Heisenberg 
model, with an XX + YY + ZZ interaction for each edge, and the antiferromagnetic XY model, with an 
XX + YY interaction for each edge. The complexities of the ground energy problems for these models are 
unknown, although some variants have been studied [3, 4, 8 ]. Recent work has established QMA-completeness 
for the antiferromagnetic XY model with coefficients that vary throughout the graph and depend on system 
size [7]. This result holds even when the graph is a triangular lattice in two dimensions. Our result is 
incomparable because while we consider general (simple) graphs, we restrict the coefficient for every edge to 
be 1. In addition, the restriction to fixed magnetization is not present in reference [ ]. 

Our main result concerns the antiferromagnetic XY model on a graph. For a given simple graph G with 
vertex set V(G) and edge set E(G ), the Hamiltonian has the form 

(1) \ X {XiXi + Y&i) 

where A, Y, Z denote Pauli matrices and a subscript indicates which qubit is acted on. Note that this 
Hamiltonian commutes with the total magnetization operator M z = X);eu(G) ^i- so it decomposes into 
sectors for each eigenvalue of M z . We prove that approximating the ground energy of this Hamiltonian in a 
sector with fixed magnetization is a QMA-complete problem. 

Our result is a natural extension of our previous work [; ]. The difference is that we previously considered 
Hamiltonians defined by graphs that may have self-loops. In that context we established QMA-completeness 
of the ground energy problem for the Bose-Hubbard model (at fixed particle number), a system of bosons 
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hopping on a graph with a repulsive on-site interaction. Then, using the relationship between hard-core 
bosons and spins, we also established QMA-completeness of a ground energy problem related to the XY 
model, where the Hamiltonian has an XX + YY term associated with each edge (as in (1)) as well as a local 
magnetic field associated with each self-loop. In this paper we present a stronger result since we consider 
only simple graphs and Hamiltonians of the form (1). 

Our proof relies heavily on machinery and results from reference [3]. We prove QMA-hardness of the 
ground energy problem for the Bose-Hubbard model on a simple graph at fixed particle number. Our 
starting point is the previous QMA-completeness result [3] which pertains to graphs G with self-loops. 
Counterintuitively, we begin by increasing the number of self-loops: we modify a graph G from the previous 
construction to obtain a new graph G SL in which every vertex has a self-loop. The new graph G SL is formed 
by taking two copies of G, adding self-loops to each of them, and then adding some edges between the two 
copies. The modification is performed in such a way that the ground spaces of the Bose-Hubbard model on 
G and G SL (in the sector with a given number of particles) are simply related. Once we have a self-loop 
at every vertex, we then remove all self-loops to obtain another graph G' NSL with no self-loops, which was 
our goal. This removal is equivalent to subtracting a term in the Hamiltonian that is proportional to the 
identity. We thus obtain a graph G NSL with no self-loops such that approximating the ground energy of the 
Bose-Hubbard model on G NSL is as hard as approximating the ground energy of G, which is QMA-hard by 
our previous result. 

Finally, we use a reduction presented in [3] showing that an instance of this ground energy problem for 
the Bose-Hubbard model on a graph is equivalent to an instance of the ground energy problem for the XY 
model at fixed magnetization on the same graph. Since this reduction preserves the graph, we establish 
QMA-completeness for the XY model as discussed above. 

The remainder of this paper is organized as follows. In Section 2, we provide basic definitions and tools 
used in this paper. In Section 3, we review results from reference [3]. In Section 4, we describe a procedure 
that modifies graphs from the previous QMA-completeness result so that every vertex has a self-loop. We 
show a relationship between the ground energies before and after the modification. Finally, we remove all the 
self-loops from the modified graph, giving only a constant overall energy shift in the associated Hamiltonian. 
In Section 5, we use this strategy to establish that the ground energy problem for the antiferromagnetic XY 
model (on simple graphs, at fixed magnetization) is QMA-complete. 

2. Preliminaries 

In this paper, G denotes a graph with vertex set V(G), edge set E(G), and adjacency matrix A(G). Later 
we will be interested in the case where G is a simple graph, but for now we allow the possibility that it has 
at most one self-loop per vertex. In other words, A(G ) can be any symmetric 0-1 matrix. 

2.1. The antiferromagnetic XY model on a graph. We define the antiferromagnetic XY model with 
local magnetic fields on G to be the |H(G)|-qubit, two-local Hamiltonian 


(2) 

o G = £ 

(|01)<10| + |10)<01|)u + 

E i 1 )^ 


A(G)y= 1, 


A{G ) U =1 

(3) 

= E 

XiXj + YiYj \ - 

2 + E 

1 -Zi 

2 ‘ 


A(G)y=l, 

A(G) U = 1 



Note that the second term is only present if the graph G has self loops; this term vanishes for simple graphs, 
giving the usual antiferromagnetic XY model. It is easy to see that this Hamiltonian (either with or without 
the second, local magnetic field term) conserves Hamming weight. Let 

Wtjv = spamj» : 2 e {0,l} |y(G)l , wt(s) = N} 

be the subspace with Hamming weight N. We write 0n(G) for the smallest eigenvalue of Oq within the 
sector with Hamming weight N (i.e., the smallest eigenvalue of the restriction Odwtiv)- 


2.2. The Bose-Hubbard model on a graph. We now review the Hamiltonian of the Bose-Hubbard model 
on a graph, as defined in [3]. We present only the “first quantized” formulation of this model that we use in 
this paper; see [ ] for a broader discussion. 
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The Hilbert space of N distinguishable particles that live on the vertices of G is 

(C |v <. G )I) (8JV = span{|*ii 2 ... *at) : ij £ V[G)}. 

Here each register represents the location of a particle. The Hilbert space of N indistinguishable bosons on 
G is the subspace of symmetric states 


Z n (G) = span{Sym(|*i *2 ■ ■ -i N )) ■ ij £ V(G)} 


where 

Sym(|ii ? 2 ■ ■ ■ in }) = 

and Sn is the symmetric group on N elements. 
For any operator M £ write 


1 

Vm 


^ ' l*7r(l)*7r(2) • ■ ■ hr(IV)) 
TTESn 


M (i) = 1®• • •®1 ®M ® 1 ® • • • ® l . 

i— 1 N—i 


for the operator that acts on ( '<C\ V ( G '>\)® N as M on the ith register and as the identity on all other registers. 
Define 

N 

(4) H% = Y.A{G) {i) + E 

i=i keV(G) 


where 


N 

nk = E i fc )( fc i w 

i=l 


is an operator that counts the number of particles at vertex k. The Hamiltonian (4) acts on the distinguishable- 
particle Hilbert space. Since it is symmetric under permutation of the N registers, the bosonic space Zm{G) 
is an invariant subspace for Hq . The iV-particle Bose-Hubbard model is the restriction of this Hamiltonian 
to the bosonic subspace 


II, 


N 


M n \ 

I Z N (GY 


Looking at equation (4), we see that the first term (the movement, or hopping, term) has smallest eigenvalue 
equal to N times the smallest eigenvalue of A(G ), which we denote /i(G), while the second term (the 
interaction term) is positive semidefinite. Hence the smallest eigenvalue of Hq is at least that of Hq , which 
is at least Nn(G). 

It is convenient to subtract this constant to make the Hamiltonian positive semidefinite. Define 


H{G, N) = Hq — Nn(G), 


and write A ]^{G) > 0 for its smallest eigenvalue. When X 1 N (G) = 0 the 7V-particle ground space minimizes 
the energy of both terms (movement and interaction) separately, and we say the Hamiltonian is frustration 
free. 

There is a connection between the Bose-Hubbard model on a graph and the XY model on the same graph. 
Consider the subspace of bosonic iV-particle states that have zero energy for the second (interaction) term in 
(4). (For example, any frustration-free state lives in this subspace.) States in this subspace have no support 
on basis states where more than one particle occupies any vertex of the graph. This is the subspace of 
hard-core bosons. Since every vertex can be occupied by at most one particle, this subspace can be identified 
with the Hamming weight N subspace of |H(G)| qubits (each qubit represents a vertex, with basis states 
10), 11) representing unoccupied and occupied states, respectively). Thus the action of the Bose-Hubbard 
model in the subspace of hard-core bosons is equivalent to a spin model. In fact, as discussed in [3], the 
restriction of (4) to the subspace of hard-core bosons is exactly equal to the restriction of Oq (given by (3)) 
to the Hamming weight N space Wt^. 


2.3. Gate graphs. The QMA-completeness construction from reference [3] uses a class of graphs called gate 
graphs that we now review. 
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Figure 1. Diagram elements. 


2.3.1. The graph g q. A 128-vertex simple graph denoted go plays a central role in the construction. The 
graph go is defined explicitly in reference [3] by specifying its adjacency matrix A(go); here we only describe 
the properties that we use in this paper. The vertices of go are labeled by tuples 


(z,t,j) :ze{0, l}, t e [8], j e{ 0,... ,7} 
where [n] = {1,..., n}. The adjacency matrix A(go) acts on the Hilbert space 


span{ l^r) |£) |j) :z6{0,1}, t £ [8], j £ {0,..., 7}}. 

The smallest eigenvalue of A(go) is ei = — 1 — 3\/2. The corresponding eigenspace has an orthonormal basis 
given by the four states 

(5) \4>z,o) = ^(|*>(|1> + 13) + |5) + |7)) + H\z)( |2) + |8)) + HT\z)(\4) + |6»)| W > 

(6) |^,i) = \4>z,o)* 


where z £ {0,1}, 





0 


and 

h = ^X>-Da 

4=0 

Note that, depending on the value of t in the second register, the first register of the state \ip z ,o) contains 
the output of a single-qubit computation where either the identity, Hadamard, or HT gate is applied to 
the state | z). (The state of the third register is in a product state with the first two, and is somewhat 
uninteresting; this register exists for technical reasons.) 


2.3.2. Gate graphs and gate diagrams. A gate diagram is a schematic representation of a gate graph. To 
define gate graphs, we first define gate diagrams and then describe how a graph is associated with each of 
them. 

The simplest gate diagrams are shown in Figure 1. These three basic gate diagrams are also called diagram 
elements. The diagram elements each represent the same gate graph which is just the 128-vertex graph go 
described above. Each diagram element has a unitary label which is either 1, H, or HT , as well as a set 
of eight circles that we call nodes. A node labeled (z, t) is associated with the eight vertices of go labeled 
(*, t,j) with j £ {0,..., 7}. Note that only half of the 16 possible nodes (z, t) appear in a given diagram 
element. Moreover, the half that does appear depends on the unitary label 1, H, or HT. 

A gate diagram is constructed by taking a set of diagram elements and adding edges between some pairs 
of nodes and self-loops to other nodes. Each node may have a self loop or an incident edge but never both 
(and never more than one edge or self-loop). If the gate diagram has R diagram elements, then each node 
can be labeled (q, z, t) with q £ [R], z £ {0,1} and t £ [8]. We write S for the set of nodes (q, z, t ) that have 
self-loops attached and we write £ for the set of pairs of nodes {(g, z, t), (q', z’. t 1 )} that are connected by an 
edge. 

A gate diagram with R diagram elements, self-loop set S , and edge set £ is associated with a gate graph 
G as follows. The vertex set of G is just R copies of the vertex set of go, with vertices labeled (g, z, t,j ) with 
q £ [R], z £ {0,1}, t £ [8], and j £ {0,..., 7}. The adjacency matrix A(G) acts on the Hilbert space 

span{|g)|£)|t)|j) : q £ [R], z £ {0,1}, t £ [8], j £ {0,... ,7}} 
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and is defined by 

(7) A{G) = lqr ® A(go) + hs + h £ 

( 8 ) h s = \<l, z ,t){<l, z A ® 1 j 

(q,z,t)£S 

(9) h £ = ^ (W,z,t)+ \q',z',t'))((q,z,t\ +{q',z',t'\)0lj 

{(? ,z,t),(q',z',t')}ee 

where we write l g and lj for the identity operators on the first and fourth registers, respectively. Note that 
hs and h £ are always positive semidefinite, so the smallest eigenvalue of A(G) is lower bounded by ei, the 
smallest eigenvalue of A(g o). When these quantities are equal, we say G is an ei-gate graph. 

Definition 1. A gate graph G is an ei-gate graph if its smallest eigenvalue is e\ = — 1 — 3y/2. 

Suppose |x) is a ground state of the adjacency matrix of an ei-gate graph G. For each q £ [i?], define 

Wz,a) = \Q)\i>z,a) 

(where the states | ip z ,a) are defined in (5)-(6)) and note (from equation (7)) that |x) G (V, where 

( 10 ) y = span{|V>|, a ) :q€[R], z, a € { 0 , 1 }} 
and h s |x) = h £ |x) = 0 . 


2.4. Spectral bounds for positive semidefinite matrices. Throughout our proof, we bound eigenvalue 
gaps of positive semidefinite matrices. For a positive semidefinite matrix A, let "/(A) denote the smallest 
nonzero eigenvalue. 

Fact 1. Let Ha and Hb be positive semidefinite matrices. Let Ha have nonempty nullspace S. Then 

l{H A + H B )< 1 (H B \s). 

Proof. Let \rf) be an eigenvector of H B \s with eigenvalue "/(H B \s). 

Suppose first that the nullspace of Ha + H B is nonempty. In this case the nullspace of Ha + H B is equal 
to the nullspace of H B \$ (since both Ha and H B are positive semidefinite), and hence 1 77 ) is orthogonal to 
this space. Thus 

7 (Ha + H b ) < (r?| H A + H B \rj) = (rj\ H B \tj) = ^/{H B \ S ). 

If instead the nullspace of H A + H B is empty, then 7 (H A + H B ) is the smallest eigenvalue of H A + H B 
and again is variationally upper bounded by (77I H A + H B \if) = (r)\ H b \rj) = 7 (/ 7 r|s). □ 


A version of the following lemma was used but not explicitly stated in reference [(>]. We gave a proof of 
this “Nullspace Projection Lemma” in [•'!] and used it extensively in that work. We recently became aware 
of another work [ ] that proves a slightly stronger bound than the one from [3]. We quote the better bound 
here. 


Lemma 1 (Nullspace Projection Lemma [1]). Let Ha,H b > 0. Suppose the nullspace S of Ha is 
nonempty and 

l(H B \s) > c > 0 and "/{Ha) > d > 0 . 


Then 


( 11 ) 


7 {Ha + H b ) > 


cd 

d+WHsW 


3. Previous results 

In this Section we summarize results from reference [.' ] that are used in this paper. 

In [ ] we considered the problem of approximating the ground energy of the Bose-Hubbard model on 
a graph at fixed particle number. We also considered a special case of this problem where the goal is to 
determine whether the Hamiltonian is close to being frustration free (recall that in our setting, frustration 
freeness means A( V (G) = 0). 
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Problem 1 (a- Frustration-Free Bose-Hubbard Hamiltonian). We are given a K -vertex graph 
G, a number of particles N < K, and a precision parameter e = y. The integer T > 4 K is provided in 
unary; the graph is specified by its adjacency matrix, which can be any K x K symmetric 0-1 matrix. 
We are promised that either A] V (G) < e“ (yes instance) or A) V (G) > e + e“ (no instance) and we are 
asked to decide which is the case. 


Here T is given in unary so that the input size scales linearly with T. An algorithm for this problem is 
considered efficient if it uses resources polynomial in i = T. 

The positive integer a in the above definition parameterizes how close to frustration free the Hamiltonian 
is in the yes case. This definition slightly generalizes one presented in [3]; the definition of the “Frustration- 
Free Bose-Hubbard Hamiltonian” problem given in that paper corresponds to the choice a = 3. Here it 
is convenient to explicitly define computational problems for each positive integer a. The proof of QMA- 
completeness presented in reference [3] applies to each a (as noted on page 11 of [3]). 

Theorem 1 ([3]). For any positive integer a, the problem a-Frustration-Free Bose-Hubbard Hamiltonian is 
QMA-complete. 

In fact, the results of reference [3] can be used to show that each of these problems is QMA-hard when 
restricted to a certain subset of instances. In particular, 

Corollary 1. For any positive integer a, the problem a-Frustration-free Bose-Hubbard Hamiltonian remains 
QMA-hard if we additionally promise that the graph G is an ei -gate graph described by a gate diagram with 
R diagram elements, satisfying 

( 12 ) l(A(G)- ei )>^ 

where Go is an absolute constant. 

In Appendix A we prove this Corollary using the results of [3]. 

Using a reduction based on the connection between spins and hard-core bosons, we established QMA- 
completeness of the following problem. 


Problem 2 (XY Hamiltonian with local magnetic fields). We are given a A-vertex graph G, an 
integer N < K, a real number c, and a precision parameter e = y. The positive integer T is provided in 
unary; the graph is specified by its adjacency matrix, which can be any K x K symmetric 0-1 matrix. 
We are promised that either 9n (G) < c (yes instance) or else On (G) > c + e (no instance) and we are 
asked to decide which is the case. 


Theorem 2 ([3]). XY Hamiltonian with local magnetic fields is QMA-complete. Moreover, there is a direct 
reduction that maps an instance of 3 -Frustration-Free Bose-Hubbard Hamiltonian specified by G, N, and e, 
to an instance of XY Hamiltonian specified by G, N, e' = j, and c = Np{G) + with the same solution. 

This reduction is presented in Appendix B of [3]. For our purposes it is crucial that the reduction does not 
change the underlying interaction graph. In this paper we show that the a-Frustration-Free Bose-Hubbard 
Hamiltonian on simple graphs is QMA-complete, and then we use the above reduction to show that the XY 
model (on simple graphs, i.e., without local magnetic fields) is QMA-complete. 

4. Adding Self-Loops 

In general, a gate graph is not a simple graph since it may have self-loops. From equations (8) and (9), 
we see that self-loops in the gate graph G arise from both self-loops and edges in its gate diagram. An edge 
or a self-loop in the gate diagram is associated with 16 self-loops or 8 self-loops in G, respectively. In this 
Section we describe a mapping from any ei-gate graph G to a modified graph G SL . The graph G' SL is not a 
gate graph. It is designed so that it has a self-loop on each of its vertices. We prove that certain properties of 
G sl are related to those of G. Ultimately our goal is to establish a relationship between the ground energies 
of the Bose-Hubbard models on these graphs. 
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4.1. Definition of G SL . Consider an ei-gate graph G described by a gate diagram with R diagram elements 
and edge and self-loop sets £ and S , respectively. 

Define A f to be the set of vertices without self-loops in G, i.e., 

(13) Af = {(q,z,t, j) : j £ {0,... ,7}, (q,z,t) £ S , and \/(q', z',t'), {( q,z,t ), <£ £}. 

Note that Af contains (q, z, t , j) if and only if it contains ( q , z, t, i) for all i £ [8], and further that for each q, 
there exists some t* such that {q, z , t*,j) £ Af for all 2 £ {0,1} and j £ {0,..., 7}. 

The vertex set of G SL is two copies of the vertex set of G. We label the vertices of G SL as 

V(G sl ) = {(q,z,t,j,d) : q £ [.R], t £ [8], j £ {0 ,... ,7}, z,d £ {0,1}}. 

We define G SL by its adjacency matrix 

(14) A(G SL ) =^(G)®i d + 2n A 7®n + 
where 

IW = E l«>H n + = l+X+l = |(|0) +11))«0| + <i|). 

The first term of (14) is just two copies of G; the second term adds edges between the two copies as well as 
self-loops to both copies. Note that every vertex of G SL has a self-loop. 

4.2. Relationship between the adjacency matrices of G and G SL . Since A(G SL ) commutes with 
lq 0 l z 0 It <8 1 j 0 n+, there is an eigenbasis for A(G SL ) where each vector is of the form \<j>) |+) or \(f>) |—). 
For states of the latter form, the second term in (14) vanishes. From this we see that, if \if) is in the ground 
space of A(G ), then | ijj) | —) is in the ground space of A(G SL ) and its ground energy is e\. We now prove 
that these states are a basis for the ground space. 

Lemma 2. Let G be an ei-gate graph and let J- be the ground space of A{G). Let 

(15) T± = ||a)|±) : |a) £ T}. 

Then the ground space of A(G SL ) is J--, and furthermore, 

(16) (2IIy ®n + )|^. + > i. 

Proof. To prove (15) it suffices to show that no state of the form \/3) |+) is in the ground space of A(G SL ) 
(since A(G SL ) commutes with l q 0 l z 0 l t 0 1 j 0 II_|_). 

Suppose (to reach a contradiction) that |/3) |+) is in the ground space. Since the ground energy of A(G SL ) 
is equal to ei, we must have 

(P\(+\A(G^)\/3)\+) = (f3\A(G)\/3) = e 1 

and hence 

(17) (P\ntf\P) = 0. 

Now, since |/3) is in the ei-energy ground space of A(G ), it is contained in the space y defined in equation 
(10). We now show that 

(is) n*\ y > 

which contradicts (17); this will show that no such state \P) |+) exists. It also establishes (16). 

To establish (18), first observe that for q s we have (if^. b \TVjq-\ipl a ) = 0, so the operator Iljv/jy is block 
diagonal with a block for each q £ [i?] (each block is a 4 x 4 principal submatrix corresponding to the 
subspace spanned by the states \if1 a ) with z,a £ {0,1}). It is therefore sufficient to establish (18) for each 
block individually. 

Focus on the block labeled by some qo £ [i?]. Now we use the fact about Af that is noted in the text 
following equation (13), namely, that there exists some t* such that (qo,z,t*,j) £ Af for all 2 £ {0,1} and 
j £ {0,..., 7}. Using this fact we can write 

Ily = |9o)(go| ® lz 0 |t*)(i*| ®1 j+B 
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for some positive semidefinite operator B. To finish the proof of equation (18), we show that the first term 
on the right-hand side is strictly positive within the block labeled by qo , which follows from 

«“bl(ko)(go| ® U ® |t*)<t*| ® IjOI^o) = {^x,b\ l z ® ® 1 j\lpz,a) = ^5z, x S a , b , 

where in the last equality we used equations (5) and ( 6 ). □ 

Lemma 3. Let G be an e\-gate graph. Then 

704 (G SL )-ei)>G 7 04 (G)-ei) 

where C is an absolute constant. 

Proof. We use the Nullspace Projection Lemma. Write H(G SL ) — ei = Ha + Hb as a sum of two positive 
semidefinite operators 

h a = a(g)® i-ei h b = 2ny®n + . 

Note that ||.Hb|| = 2 since ny <gt II + is a projector. We will also need a bound on 7 (i?s|s) where S is the 
nullspace of H A - Note that S = J- + + T_ where T± is defined in equation (15). By Lemma 2, the nullspace 
°f Hb\s is -+-• Using this fact, we see that 7 (i?s|s) is equal to the smallest eigenvalue of Hb within the 
space Jq_, which we bound using equation (16): 

t(H b \s) > 2' 

Now applying the Nullspace Projection Lemma, we get 
(19) 7 (A(G sl ) - ei ) > 


17 ( A (G) ® 1 - ei) 


> 


h{A(G) ® 1 - ei) 


7 (A(G) ® 1 — ei) + 2 ||A(G)|| + |ei| + 2 

To complete the proof, we show that ||Al(G)|| is upper bounded by an absolute constant. Looking at the 
general expression (7) for the adjacency matrix of a gate graph, we see that 

P(G)|| < \\A(g 0 ) || + IIMI + IIMI < P(ffo)|| + 1 + 2. 

This completes the proof: the right-hand side of this expression is an absolute constant since go is a fixed 
128-vertex graph. □ 

4.3. Relationship between the Bose-Hubbard models on G and G SL . We begin by defining a linear 
map from the Hilbert space (C++ G - ) I)® JV of N distinguishable particles on G to the corresponding space 
((q|v(G )|)0IV £ or G SL_ 'pjjjg lna p i s defined by its action on basis states as follows: 

|*l)|*2) ■ ■ ■ |*iv) ^ |*li —) 1*2, -) • ■ ■ |*IV, -) 

where |u,—} = ^g(|u,0) — |u, 1)) is a superposition of the two vertices (v, 0) and (y, 1) of G SL that are 
associated with vertex v £ V(G). 

For a state |0) £ (<C.\ V ( G )\^ N we write | 0) for its image under this mapping, i.e., 

|0) = a x\ x \)\ x 2) ■ ■-\X N ) |0) = ^2 Ux\Xl,-)\x 2 ,-) ■ ■-\X N ,-)- 

xeV(G) N xEV(G) n 


Clearly, if |0) is normalized then so is \cj>). Furthermore, if 
permutations of the N registers) then |0) £ Zjv(G SL ). 


£ Zn(G) (i.e., if it is symmetric under 


Lemma 4. Let 


£ (cl y ( G )l)® Jv . Then 

01 n w {n w - 1 ) | 0 ) = 


1 


y: n v (n v - 1) 

Kvev(G) t 


\»GV(G sl ) / 

Proof. Writing (u,0) and (v, 1) for the two vertices of U(G SL ) corresponding to a vertex v £ V(G), we have 


N 


T, h w (n w -1 ) = y n v y{h v y- 1) n v y = y\v,d)(v,d\ 
i«ev(G SL ) wev(G),de{o,i} »=i 


(0 




COMPLEXITY OF THE XY ANTIFERROMAGNET AT FIXED MAGNETIZATION 


9 


We can then rewrite each term in the sum as 

n v> d(n v> d - i) = (l u > d l 0) ) (l u > d )( v > d \ (j } ) ■ 

Using this expression (and a similar expression for h v (n v — 1)), we see that 
(<?| n v ,d{n v ,d - = ($| (\v, d)(v, d\ {j 

= ((-i^)(^i-)) 2 (^i (i w >H 0) ) (V>h° } ) w 

= -(</>|n„(h„ - 1)|^). 

Summing both sides over d G {0,1} and v € V(G) gives the claimed result. □ 

Lemma 5. Let G be an ei-gate graph. Then 

(20) A)v(G) < a => A]v(G SL )<^a. 

If in addition G is described by a gate diagram with R diagram elements and satisfies 7 (A(G) — ei) > ^ 
for some absolute constant Co, then 

(21) Ajv(G) >b>0 =► A)v(G sl ) > ^r 2 b 

where C\ is an absolute constant. 

Proof. Let | <f>) G Zn(G) be a normalized state with minimal energy for H(G,N), i.e., {(f>\H(G, N)\<f) = 
A jy(G). The normalized state \<fi) G Zn{G sl ) satisfies 

A(G sl )^|0) = (A(G) ® 1„ + 2T\ m ® n + )«|^) = (A(G) <g> l d ) (i) | J) 


for each i G [iV], so 


N 


N 


I £(a (g sl ) - ei )W|^> = w ^(A(g) ® i d - ei )W| 

2=1 2=1 

N 

= ^|^(^(G)- ei )«|^ 


2=1 
/ N 


£(>*(G)-e 1 )« + £ h v (h v — 1) 

\i=i veV(G) ) 


= a 


where in going from the second to the third line we used the fact that h v (h v — 1) is positive semidefinite. 
Using this inequality and Lemma 4, we have 

A]v(G SL )< &\H(G Sh ,NM 

<a+(cf >| ^2 n w (n w - l)\(j>) 

™ev(G SL ) 

= a +\{<t>\ n v (n v - l)\(j)) 

vGV(G) 

<a+\^\H{G,NM 

3 

“ 2 °’ 


which establishes equation (20). 



10 


COMPLEXITY OF THE XY ANTIFERROMAGNET AT FIXED MAGNETIZATION 


Now suppose A] V (G) > b > 0 and y(H(G) — ei) > and consider the second part of the Lemma. Note 
that in this case H (G, N) has no nullspace, so A ]^(G) = 7 (H (G, N)). We write H (G SL , N) = Ha + Hb with 
positive semidefinite operators 


N 


H A = (G sl ) - ei ) 


(<) 


•Hat(G sl ) 


Hb — 1 )| 2 jv ( G sl) 

™GV(G SL ) 


and we use the Nullspace Projection Lemma. 

To get a bound on 7 (Ha), first note that every eigenvalue of Ha is also an eigenvalue of the operator 


( 22 ) 


N 

J2(A(G SL ) - e^W 


i= 1 

(without the restriction to Zn(G sl )), since this operator is permutation symmetric and preserves the sym¬ 
metric subspace. Hence the smallest nonzero eigenvalue of Ha is at least that of (22) and 

/ N 

(23) 7(#a)> 7 ^(A(G sl ) - ei) (i) 

\i=i 

where G is an absolute constant (in the second step we used the fact that ® I + I ® M) = 7 (M) for 
any Hermitian matrix M with smallest eigenvalue zero, and in the third step we used Lemma 3). 

We also need a bound on 'y(Hs\s)i where S is the nullspace of Ha- Letting T be the nullspace of 


= 7 (H(G sl ) - ei ) > G 7 (H(G) - ei) > 


GG 0 

R 3 


£(A(G)- ei )« 


i=1 


Z N (G) 


and using Lemma 2, we see that S is equal to the image of T under the mapping |(/>) 1 —> \<j)). Using this fact 
and Lemma 4, we get 

Hb\s = ^ n v (hv-l)\ T . 

vev(G) 


Since H(G,N) has no nullspace, neither does the operator on the right-hand side of this equation. Hence 
Hb\s > 0 and A] V (G SL ) = 7 (iL(G SL , iV)). Furthermore 


(24) 


l(H B \s) = -jl 


T. n v (n v — 1)| 

\vG V(G) 


>^(H(G,N))> b - 


where in the first inequality we used Fact 1. 

Now using the bounds (23) and (24) along with the fact that ||Hb|| < N 2 and applying the Nullspace 
Projection Lemma, we get 


A] v (G SL )= 7 (^(G SL ,iV))> 


CCob 

2 R 3 

g § L + N 2 


> 


CCob 

~ 2 W 


(GGo + l)iv 2 ■ 


Thus the result follows with Gi = 2 cc 0 + 2 


□ 


5. Removing self-loops 

Our goal is to consider simple graphs, but so far we have described a method for mapping an ei-gate 
graph G to a graph G SL with self-loops on every vertex. We now remove all the self loops from G SL to obtain 
a simple graph G' NSL . The adjacency matrix of this graph is 

H(G nsl ) = A(G sl ) - I 

and it has smallest eigenvalue /r(G NSL ) = /k(G sl ) — 1. 

Now consider the TV-particle Bose-Hubbard model on G NSL . We have 

N 

Hg ksl = ^ H(G nsl ) w + y h k {fi k - 1) = H% 3L - N, 

»=1 fcGV(G NSL ) 
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H(G nsl ,N) = Hg NSL - iV M (G NSL ) = (Hg* -N)- iV(/i(G SL ) - l) = H(G SL ,7V). 
In particular, the smallest eigenvalues of these two Hamiltonians are equal: 

(25) A] v (G nsl ) = A] v (G sl ). 

We now use this relationship to show that the following problem is QMA-complete. 


Problem 3 (a- Frustration-Free Bose-Hubbard Hamiltonian on simple graphs). We are given 
a AT-vertex simple graph G, a number of particles N < AT, and a precision parameter e = y. The integer 
T > 4A' is provided in unary. We are promised that either A ]^(G) < e a (yes instance) or A) V (G) > e + e“ 
(no instance) and we are asked to decide which is the case. 


Theorem 3. For any positive integer a, the problem a-Frustration-Free Bose-Hubbard Hamiltonian on 
simple graphs is QMA-complete. 


Proof. Let a be a fixed positive integer. The problem is clearly contained in QMA since it is a special case of 
the QMA-complete problem a-Frustration-Free Bose-Hubbard Hamiltonian. To show that it is QMA-hard, 
we reduce from another QMA-hard problem. Let /3 = 8a and define “Problem A” to be the special case 
of /3-Frustration-Free Bose-Hubbard Hamiltonian where the graph G is promised to be an ei-gate graph 
satisfying (12). Corollary 1 implies that Problem A is QMA-hard. We provide a reduction from Problem A 
to a-Frustration-Free Bose-Hubbard Hamiltonian on simple graphs. 

Let an instance of Problem A be given, specified by G, N, and e. We assume that e is smaller than some 
absolute constant G. We show that any such instance of Problem A has the same solution as the instance 
of a-Frustration-Free Bose-Hubbard Hamiltonian on the simple graph G' NSL , with number of particles N 
and precision parameter e' = e 7 . This is sufficient to prove QMA-hardness of a-Frustration-Free Bose- 
Hubbard Hamiltonian on simple graphs since there are only finitely many instances of Problem A (and of 
/3-Frustration-Free Bose-Hubbard Hamiltonian) with e lower bounded by a constant. 

First we check that G NSL , N, and e' satisfy the conditions in the definition of the problem, i.e., that 
they specify a valid instance of a-Frustration-Free Bose-Hubbard Hamiltonian. Let K = |H(G)| and K' = 
|F(G NSL )|. Then K' = 2 A', and since e = ^ we have 



with T' = T 7 > 8K = 4AT'. Furthermore, N < K < K'. Therefore G NSL , N, and e' specify a valid instance. 

Next we show that, provided e is smaller than some constant G, the instance of a-Frustration-Free Bose- 
Hubbard Hamiltonian defined by G NSL , N, and e' has the same solution as the instance of Problem A defined 
by G, N, and e. 

Suppose first that the instance of Problem A is a yes instance. Then A]y(G) < t 13 . Applying the first part 
of Lemma 5 and using (25), we get 


A )v(G nsl 



-e 8a < e 7 “ = (e') a 
2 - v ; 


where in the second-to-last inequality we used the fact that |e“ < §e < 1, which follows from the fact that 
e -1 is an integer greater than 1. So in this case G NSL , N, and e' specify a yes instance of a-Frustration-Free 
Bose-Hubbard Hamiltonian. 

Next suppose that the instance of Problem A is a no instance, so A)y(G) > e + eP. Applying the second 
part of Lemma 5 and using (25), we get 


A V(G“ SL )>A^( e + e'>) 

where D is an absolute constant and R is the number of diagram elements in G. Noting that R < K, N < K, 
and e < we see that 

A]v(G NSL ) > D'(e 6 + e 5+/3 ) 
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where D' is another absolute constant. Now, provided 2e < D' (which holds for sufficiently large problem 
size K), we have 

A]v(G nsl ) > 2e 7 > e 7 + (e 7 )“ = e' + (e')“, 

which shows that G NSL , N, and e' specify a no instance of a-Frustration-Free Bose-Hubbard Hamiltonian. □ 
Finally, we apply this result to the XY Hamiltonian problem. 


Problem 4 (XY Hamiltonian on simple graphs). We are given a A'-vertex simple graph G, an 
integer N < K , a real number c, and a precision parameter e = A. The positive integer T is provided 
in unary. We are promised that either 9n(G ) < c (yes instance) or else 9x(G) > c + e (no instance) and 
we are asked to decide which is the case. 


QMA-hardness of this problem now follows directly from the fact that 3-Frustration-Free Bose-Hubbard 
Hamiltonian on simple graphs is QMA-hard (from Theorem 3) and the reduction from Theorem 2. 

Theorem 4. XY Hamiltonian on simple graphs is QMA-complete. 
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Appendix A. Proof of Corollary 1 

The strategy used in [3] to establish QMA-hardness is to reduce quantum circuit satisfiability to Frustration- 
Free Bose-Hubbard Hamiltonian using two steps. First, starting from a QMA verification circuit Cx for an 
instance X of a problem in QMA, it was shown how to construct an ei-gate graph Gx • (The explicit 
construction of Gx is detailed in Section 6.2 of [3], and the fact that it is an ei-gate graph is proven in 
Lemma 11 of that paper.) Second, the graph Gx (along with some extra information encoding “occupancy 
constraints”) is used to construct another ei-gate graph G^ (this construction is detailed in Appendix C of 
[3]). The class of instances of Frustration-Free Bose-Hubbard Hamiltonian that are shown to be QMA-hard 
all use graphs Gy of this form. Thus, to prove Corollary 1, it suffices to show that 

T{A{Gx)~ei) > 


(26) 
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where R is the number of diagram elements in G§- and Co is an absolute constant. 

To prove the following results, we rely on some facts established in [3] (in particular, see Sections 5, 6.2, 
7.1, and Appendix E.4). In this appendix we analyze the spectral gaps of Gx and G^ and thereby establish 
(26). 

Lemma 6 . 7 (A(Gx) — ei) > D, where D is an absolute constant. 

Proof. Section 6.2 of [3] describes in detail how the graph G\ is constructed from the verification circuit. 
Its adjacency matrix can be written 

A{G X ) = A{Gi) + h e > + h S ‘ 

where G\ is an ei-gate graph with many components, each of constant size. Here /i£< and hs' are of the 
form given in equations ( 8 ) and (9), but in this case, £' and S' are subsets all the edges and self-loops, 
respectively, in the gate diagram for Gx- In particular, the gate diagram for G\ contains some of the edges 
and self-loops in the gate diagram for Gx, and the rest are included in the sets £' and S'. 

Note that in reference [3] the terms we refer to as hs 1 and h$' above are denoted by 

n 

( 27 ) fig* + hs r — h\ + /z-2 E /lin,i T /lout- 

i=n in+1 

Each of the terms on the right-hand side corresponds to a subset of the edges or self-loops in the gate diagram 
for Gx; in reference [3] it was convenient to further subdivide the sets £' and S' , but it will not be necessary 
to review the details of this here. 

Every component of A(Gi) is one of 5 fixed graphs, each with a constant number of vertices (four of these 
5 possible graphs are called two-qubit gate gadgets and the fifth is called the boundary gadget-see Figures 
5.2-5.3 of [3]). The smallest nonzero eigenvalue satisfies 

(28) 7 (^(Gr) - er) > c, 

where c is an absolute constant equal to the smallest nonzero eigenvalue of one of these 5 graphs. 

A basis for the ei-energy ground space of A(Gi), which we denote by S in the following, is given by the 
set of states {| a ) : L £ £} defined before Lemma 11 in [3]. In this basis for S one can compute the matrix 
elements of the operator 

(29) {hs’+hs’) \ s 

These matrix elements are computed in reference [3]. Specifically (referring to the decomposition on the 
right-hand side of equation (27)), the matrix elements of h\\s are given in equations (E29)-(E34), those of 
^2 |s are given in equations (E39)-(E40), those of ft.i nj i|s appear in the equation preceding before (E45), and 
those of /lout|s are given in the second equation of Section E.4.4. 

Here we use the following fact about the operator (29), which follows from the equations of reference [3] 
listed above: it is block diagonal in the basis {| a ) : L £ C} with blocks of size at most two. Furthermore, 
each of these blocks is one of a fixed set of 2 x 2 (or lxl) matrices, and therefore the smallest nonzero 
eigenvalue of (29), equal to the smallest nonzero eigenvalue of one of these blocks, is lower bounded by a 
fixed constant that we denote c. Now we apply the Nullspace Projection Lemma using this fact, equation 
(28), and the fact that 

II hs’ + h$’ || < || hgf || + ||/i<s' || <2 + 1 . 

We get 

7 (A(G x )-e 1 )>-^~, 
c + 3 

which completes the proof. □ 

The following lemma uses a construction from Appendix C.3 of [3]. 

Lemma 7. 7 (A(Gy) — ei) > ^, where R is the number of diagram elements in the gate diagram for G ^ 
and Gq is an absolute constant. 
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Proof. The adjacency matrix of the graph can be written 

^(^x) = + hg o + hg 0 , 

where G* is an ei-gate graph (it is described in equation (C.8) of [3] and is used to enforce the occupancy 
constraints, although we will not need its explicit form) and hg 0 and hg 0 are of the form in equations (9) and 
(8), respectively, for some set of edges and self-loops in the gate diagram for G^-- States with low energy for 
the graph Gy are related to those with low energy for Gx in the following way. 

Since Gx is an ei-gate graph, its adjacency matrix can also be written as 


A(Gx) — 1 <j ® ^ 4 ( 5 o) + hg + hg. 

Let Q be the ei-energy ground space of the first term in the above equation (i.e., l q (g> A(go)), and let Q A 
be the ei-energy ground space of A(G a ). The graph G* is constructed so that these two spaces have the 
same dimension, and furthermore 

R' odd 


( hg 0 + hg 0 ) | q & = (hg + lig ) | q ■ 


3fi'+2 

3fl'-l 


R' even 


where R' is the number of diagram elements in the gate diagram for Gx (this is equation C.51 of [3]). It is 
also proven in Lemma 19 of [3] that 

- ei ) > 1 


(30) 

Using this equation, we get 


(30 R ') 2 


(31) 


l((hg 0 + hg 0 )|) > ^77 ((hg + h s )\ Q ) > ^7 (A(G X ) - ei) > ^7 


where in the second inequality we used Fact 1 and where D is the absolute constant from Lemma 6. 

We now apply the Nullspace Projection Lemma with the decomposition A(G^-) — e\= Ha + Hb, where 
Ha = A(Gx) — ei and Hb = hg 0 + hg 0 . Note that ||U_b|| < ||Ii£|| + ||frs|| < 3 . Applying the Nullspace 
Projection Lemma and using (31) and (30), we get 

" Co 


l{ A {Gx) - ei) > 


4iJ'(30iJ') 2 ^ 

+~3 - W 


1 

J30WJZ 


where Go is an absolute constant. To complete the proof, we use the fact that R (the number of diagram 
elements in G^) is greater than or equal to R' (the number of diagram elements in Gx), which is a general 
property of the mapping G G n as described in Appendix C of [3]. □ 








